Charge rearrangement and screening in a quantum point contact 
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Compressibility measurements are performed on a quantum point contact (QPC). Screening due 
to mobile charges in the QPC is measured quantitatively, using a second point contact. These 
measurements are performed from pinch-off through the opening of the first few modes in the QPC. 
While the measured signal closely matches a Thomas-Fermi-Poisson prediction, deviations from 
the classical behavior are apparent near the openings of the different modes. Density functional 
calculations attribute the deviations to a combination of a diverging density of states at the opening 
of each one-dimensional mode and exchange interaction, which is strongest for the first mode. 

PACS numbers: 73.21.Hb,73.23.-b,73.43.Fj,73.61.Ey,71.45.Gm,73.43.Cd 



The simplest of mesoscopic systems is a point contact, 
a narrow constriction between two electron reservoirs. 
Conductance measurements through such a quantum 
point contact (QPC) reveal steps in units of Go = 2e 2 /h. 
The physics of these steps is well understood: for a QPC 
adiabatically connected to the reservoirs, the transmis- 
sion coefficient of each mode is either zero or one. The 
number of such modes increases, as the width of the QPC 
is increased, leading according to the Landauer formula 
to a series of quantized steps in conductance. Such con- 
ductance measurements, however, also reveal our limited 
knowledge of the physics of QPCs at low densities. As 
a QPC is just being opened up, its conductance pauses 
around 0.7 x Go, before rising to the first full-channel 
plateau. This "0.7 structure" has been one of the prime 
puzzles in mesoscopic physics jl,]. Longer ID wires show 
a similar, and likely related, structure at 0.5 x Go 0]. 
These features have been variously attributed to sponta- 
neous spin polarization in the QPC , to Luttinger 
liquid behavior j(| , or to the formation of a localized mo- 
ment at the QPC [7j, |8|, [9(, together with the resulting 
Kondo effect. 

To shed more light on the physics of QPCs near pinch- 
off, we report here measurements of the compressibility 
of electrons in the channel of the QPC. To measure com- 
pressibility, one applies a potential to an electrode on one 
side of a structure of interest, and measures the potential 
on the other side. Compressibility measurements have 
helped elucidate how charge carriers arrange themselves 
in two dimensions when kinetic energy is dominated by 
interactions, whether in the quantum Hall regime [l(J or 
at low density and high effective mass Augmented 
by local electrostatic detectors, such measurements have 
produced striking images of how individual carriers local- 
ize in these same regimes Compressibility measure- 
ments using nanofabricated electrostatic detectors such 
as QPCs have also become a standard tool to probe 
transitions between charge states of a quantum dot or 
two coupled quantum dots 3, 14, In this letter we 



present analogous measurements, using a QPC as an elec- 
trostatic detector but measuring in this case the charge 
configuration of a second QPC instead of a quantum dot, 
c.f. Ref. Q. 

Compressibility of perfect ID systems is a textbook 
problem. The diverging density of states should give 
strong compressibility at low carrier density. In a multi- 
mode ID system, one would expect that a similar signa- 
ture should occur at the opening of each mode. In con- 
trast, both experiment and numerical simulation of our 
short ID wire show that the enhancement in compress- 
ibility at the opening of the first mode is significantly 
different from that at the opening of higher modes. In 
retrospect, this might have been expected because at the 
opening of the first mode the total density is small and 
thus exchange effects are also important. In addition, 
the distinctive "0.7" transport feature might have had 
a counterpart in compressibility. As it turns out, the 
simulated compressibility feature associated with the 0.7 
regime was too weak to observe experimentally with our 
current sensitivity. 

The devices in this experiment were fabricated on 
a GaAs/AlGaAs heterostructure, containing a two- 
dimensional electron gas (2DEG) 70 nm below the sur- 
face. The electron density n s — 2 • 10 11 cm -2 and mobil- 
ity [i = 2.3 • 10 6 cm 2 /(Vs) at 4.2 K. A schematic of the 
measured devices is shown in the inset of Fig. 1(a). Two 
QPCs, with lithographic widths of 330 nm and 350 nm, 
respectively, are separated by an 80 nm wide gate. Al- 
though the device is symmetric in design, the two QPCs 
play fundamentally different roles in our experiment: the 
lefthand QPC serves as a detector, sensitive to charge re- 
arrangements in the righthand QPC. For clarity, we will 
henceforth refer to the righthand QPC as "the QPC" and 
the lefthand QPC as "the detector." Our data support 
the assumption that the primary interaction between the 
two QPCs is electrostatic. Two nominally-identical de- 
vices, each containing a QPC and a detector, were mea- 
sured in a 3 He cryostat with 300 mK base tempera- 
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ture. The extensive measurements taken on one QPC- 
detector pair are presented in this paper; consistent be- 
havior was also found after thermal cycling to room tem- 
perature. Measurements of another nominally-identical 
QPC-detector pair were used to confirm qualitative fea- 
tures of the data presented here. 

Fig. 1(a) and (b) show the linear and non- linear differ- 
ential conductance, respectively, of a QPC with the clas- 
sic signatures of 0.7 structure. The charge detector signal 
(Fig. 1(c)) is measured as the QPC is opened from pinch- 
off through the third plateau. The full range of V qpc is 
broken into 14 shorter measurements and the detector 
readjusted for each, to keep the detector conductance in 
its sensitive, near-linear regime. Each of the 14 traces 
in Fig. 1(c) covers a 100 mV range in Vq pc , with consec- 
utive traces overlapping by 50 mV. To reduce noise we 
averaged each trace over 40 measurements. No striking 
features are visible in these raw data. 

To analyze the data more carefully, we perform two 
transformations. First, we assign to each lq pc an effec- 
tive voltage, V e ff, that would produce the same change in 
the detector conductance if it were applied to the detector 
gate. This mapping is based on the measured response of 
the detector Gdet(Vdet)) Fig. 1(c) inset. The result of this 
transformation is shown in Fig. 1(c) (right axis). Second, 
the derivative D = dV c s/dV qP c eliminates the offsets be- 
tween consecutive traces. This quantity describes the 
relative coupling of the detector to y qpc and Vdct, and 
thus provides a quantitative measure of the screening of 
V qpc by mobile charges in and around the QPC. 

The derivatives of the 14 measurements shown in 
Fig. 1(c) are plotted in Fig. 1(d). The overlapping ranges 
of each curve agree, allowing us to extract a continu- 
ous D(V qpc ). The QPC conductance is superimposed 
on the measurement of D in Fig. 1(d). Three impor- 
tant features are observed. First, D increases steadily as 
V qpc becomes more negative. Second, the slope of D be- 
comes larger immediately beyond pinch-off of the QPC 
(Vqpc < —0.95 V). Third, each steep rise in the QPC con- 
ductance is accompanied by a slight dip in D, indicated 
with dashed lines in Fig. 1(d). 

These features may be qualitatively understood as fol- 
lows. The first and second features are caused by the 
reduction in screening of V qpc by the 2DEG, as Vqp C is 
made more negative and pushes the 2DEG away from the 
QPC. This recession becomes more rapid after pinch-off. 
The explanation of the dips associated with the open- 
ing of each mode is less straightforward. Intuitively, the 
divergence of the ID density of states (1DOS) of the elec- 
trons in the QPC at the entry of each additional subband 
should be accompanied by an enhancement of the screen- 
ing and a dip in D. Less obviously, enhanced screening at 
pinch-off could also result from the exchange interaction 
between conduction electrons. 

To understand the features in D both qualitatively and 
quantitatively, we first simulate the device numerically, 




Figure 1: Linear (a) and non-linear (b) conductance of the 
device sketched in inset of (a). Plateaux fall below Go due 
to finite series resistance, not subtracted off in figures. Char- 
acteristic signatures !)| of 0.7 structure are observed: high- 
bias plateau near 0.8(2e 2 /^) (dashed ellipse) and zero-bias 
anomaly (solid ellipse), c) Inset: Conductance of the detector 
QPC, used to derive the mapping Viff(Vqpc)- Optimal work- 
ing range of the detector indicated by dotted ellipse. Main 
panel: 14 charge sensing measurements covering wide range 
of Vqpc. Left axis: detector conductance Gdet- Right axis: 
V e f{ mapped from V qpc as described in text, d) Conductance 
trace of the QPC (right axis) and derivative D = dV c s / dV qpc 
of the detector data (left axis) for identical gate settings. D is 
a measure of screening of electric fields by the QPC. Dashed 
lines mark steep rises in the conductance of the QPC, coincid- 
ing with suppressions in D which indicate enhanced screening. 



using M. Stopa's SETE code [3]. The simulation cal- 
culates self-consistently the effective potential and the 
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density in the QPC as a function of voltages on the elec- 
trostatic top gates, but does not include quantum cor- 
rections associated with the ID constriction in the QPC 
[Hi]. Inputs to the simulation include the 2DEG growth 
parameters, geometry and voltages of the three gates, 
and temperature. The potential landscape of the device 
is calculated for a range of settings of V^ pc and Vdet- 

In Fig. 2(a) we plot D SETE extracted from the simu- 
lated data using a procedure analogous to that described 
for D above. The quantitative match between the sim- 
ulated £) SETE and the measured D is striking. In the 
simulation, as in the measurement, the kink in £) SETE 
occurs exactly at pinch-off [13] ■ This confirms the asso- 
ciation of the change in slope with an electrostatic effect 
of emptying the point contact's saddle potential. 

Both the measured and the simulated derivatives are 
linear above and below the kink at pinch-off. As seen in 
Fig. 2(a), the simulation is closely approximated by two 
lines intersecting at the pinch-off voltage. Fig. 2(b) has 
analogous lines to approximate the data above and below 
pinch-off and hence to emphasize the fact that the data 
drop below these guide lines at the transitions between 
conductance plateaux. The absence of such dips in the 
simulation, which includes only the classical electrostatic 
effect of the gates on the 2DEG, supports the argument 
that the modulations are caused by quantum mechanical 
effects such as a nontrivial 1DOS or exchange interaction. 

The data in Fig. 2(b) show a pronounced dip at pinch- 
off. Understanding the source of this dip and the smaller 
dips associated with opening of successive subbands re- 
quires a quantum mechanical calculation. Since such a 
calculation of the full three dimensional system is very de- 
manding, we perform a density functional theory (DFT) 
simulation of a device in which the QPC is modeled by a 
constriction in a long quantum wire that is wide enough 
to carry four spin-degenerate modes [2~ll. |22T| . In the inset 
of Fig. 2(c), the quantum wire is defined by gates marked 
V c and the constriction by gates marked V qpc . We fur- 
ther modify the dimensions and electron density to make 
the computation feasible. The width and the length of 
the simulated QPC are 250 nm and 200 nm, respectively, 
and the 2DEG is 70 nm below the surface. The donor 
density is 10 11 cm -2 . The screened potential Vdot2 is de- 
tected 200 nm beneath the center of the QPC. We simu- 
late the device for a range of Vq pc such that the QPC has 
from zero to three open subbands - the thin solid line in 
Fig. 2(c) shows the conductance. Screening of V^pc by the 
QPC affects the value of £> DFT = dV dct 2/dV qpc . We use 
the local density approximation (LDA) for the exchange- 
correlation energy of the electrons in the 2DEG, us- 
ing the known effective mass and dielectric constant for 
GaAs. Because of the modified geometry and the dif- 
ferent electron density the DFT simulation cannot re- 
produce our measured D quantitatively. However the 
simulation clearly shows a dip whenever a new subband 
opens. In agreement with the measured data, the dip 



0.12 - 



0.1 



0.06 - 



0.04 ■ 



0.12 



0.10 



0.08 



0.06 



(a) 



pinch-off 



0.6 -0.4 




0.08 



0.07 



Q 

Q 



0.06 



0.05 



-1 -0.9 


-0 


8 -0.7 
Vqpc 


-0.6 -0.5 -0.4 

(V) 


* 1 






1 1 

(c) 








v c 




1 

T qpc 1 


V 

»qpc 










v c 




LDA 








LDA, 

m* — > 2m* / 
Hartree / 




1 1 



-0.25 



-0.2 



Vqpc(V) 



-0.15 



-0. 1 



Figure 2: a) D SETE = dV c g /dV qpc extracted from numeri- 
cal simulation of the device. The red dashed lines are guides 
to the eye, showing linearity of D SETE around pinch-off. No 
modulation can be observed in the simulation, which does 
not account for the one-dimensional nature of the QPC (l9| . 

b) The shaded area indicates a pronounced dip of D around 
pinch-off. No such feature is observed in the simulated D SETE . 

c) D DFT , as calculated from DFT simulation of a modified 
QPC (inset) with a detector 200 nm beneath the plane of the 
QPC. Thick solid line: LDA with actual GaAs effective mass; 
dashed line: LDA with doubled mass; dotted line: Hartree ap- 
proximation with actual mass. The thin solid line shows the 
conductance of the QPC in the LDA approximation (right- 
hand scale). The dashed red line is again a guide to the eye, 
showing the change in slope of the detector signal at pinch- 
off within all these simulations. The simulations dip below 
this guide line at the opening of each new mode. Within the 
Hartree approximation the dips are all roughly equal in size, 
but in LDA the dip at pinch-off is substantially larger than 
the others, as observed in experiment. 
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Figure 3: Curves from Fig. 2(c) expanded and shifted verti- 
cally to highlight subtle differences among the approximations 
at the opening of the first and second subbands. The LSDA 
curve with an additional dip associated with spin rearrange- 
ment is also shown, unshifted with respect to the LDA curve. 



at the opening of the first subband is more pronounced 
than those for the second and third subbands. To esti- 
mate the importance of exchange and correlation we elim- 
inate them by making the Hartree approximation (dotted 
curve) [23j]. The dips are much weaker in this approx- 
imation (they disappear entirely in the Thomas-Fermi 
approximation) and all dips are of similar size, demon- 
strating that the larger dip at pinch-off in the LDA calcu- 
lation is dominated by exchange-correlation effects. This 
is not surprising, since the exchange-correlation contribu- 
tion to total energy at low carrier density is larger than 
the kinetic contribution: for quasi-lD systems, these two 
contributions go as n 3 / 2 and n 3 , respectively, where n is 
the linear density of electrons. The prominence of the 
dip at pinch-off relative to those for higher modes in the 
experiment (Fig. 2(b)) is thus evidence of the importance 
of exchange-correlation effects to charge distribution in a 
QPC. For higher modes, kinetic energy plays an impor- 
tant role: lowering kinetic energy by doubling effective 
mass (dashed curve, Fig. 2(c)) reduces dip area by 40%. 

To see whether spin rearrangement within the QPC 
should be detectable by our compressibility measure- 
ments, we have also performed a DFT simulation using 
the local spin density approximations (LSDA) in which 
the density of electrons with spin up is not restricted to be 
the same as that for spin down. In this case the simulated 
detector signal has an additional small dip (Fig. 3, dot- 
dashed), marking the subtle rearrangement of electron 
density associated with the formation of a quasibound 
state at the center of the QPC. This smaller dip would 
be a signature of the formation of a magnetic moment 
responsible for the Kondo effect. Only one broad dip is 
discernible in the experimental data, precluding a clear 
experimental statement about spin rearrangement. 

In conclusion, we have used a charge-sensitive detec- 
tor to measure the charge rearrangement in a QPC as 



it passes through pinch-off. The broad features in the 
detector signal closely match the predictions of a classi- 
cal electrostatic simulation of the device. However, the 
charge redistribution at the entry of ID subbands creates 
a series of dips in the detector signal. DFT calculations 
provide a good qualitative match to measurement, and 
indicate that the first dip is dominated by exchange inter- 
action between electrons. The other dips are associated 
with the divergence of the density of states at the opening 
of each subband. Similar DFT calculations which further 
allow for the formation of local (spin-degenerate) mag- 
netic polarization serve as the basis for the local moment 
Kondo scenario for 0.7 structure 0. Our measurements 
are consistent with both LDA and these LSDA calcula- 
tions, since in the calculations the spin rearrangement is 
accompanied by only a slight charge rearrangement. 
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aware of related work on carbon nanotubes 24] 
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Supplementary Information: Results of 
the DFT calculation with the potential 
detected in the 2DEG plane 



In the DFT calculation we unfortunately cannot sim- 
ulate a device in which both the measured QPC and the 
detector QPC are present. In order to perform such a 
calculation we would have to put both QPCs inside a 
very wide quantum wire. The number of modes such a 
wire would have to support would be too large for the cal- 
culation to be feasible. In the single QPC geometry we 
used in our DFT calculation we could have measured the 
potential in the 2DEG plane next to the QPC instead 
of 200 nm below the 2DEG plane as we did. It turns 
out that the features observed in the experiment and in 
our DFT calculation can also be observed if we simulate 
measuring the potential in the 2DEG plane (see Figure 
SI.) In Figure S2, the derivative D was calculated from 
the potential measured in the 2DEG plane 300 nm from 
the center of the QPC. As in the calculation presented 
in the paper, there is a large dip at pinch-off and a series 
of weaker dips corresponding to higher transverse modes 
being occupied in the QPC. Observation of these features 
in two different geometries supports our contention that 
they are qualitatively linked to the features we observe in 
the experiment. However, because the experimental ge- 
ometry is different from that of the two calculations we 
cannot claim that the calculations should quantitatively 
account for the magnitudes of the dips in the real device. 
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Figure S2: The derivative D DFT calculated from the poten- 
tial off to the side in the 2DEG plane (detector location is 
shown as a blue dot in Fig. Sib). The curves for three dif- 
ferent approximations are shown: DFT within local density 
approximation (red), DFT within local spin-density approx- 
imation (green) and the Hartree approximation (blue). The 
curves show very similar behavior to that of the calculation 
presented in the paper (Fig. 2c and Fig. 3), showing that 
the detailed geometry is not important in determining the 
qualitative features of the compressibility measurement, as 
long as the system under study lies between the gate that 
is varied and the detector. The mag nitude of D DFT , how- 
ever, is smaller here because the right gate electrode effec- 
tively screens the potential. In the experiment, the electrode 
separating the two QPCs is very thin and is thus less effective 
in screening the potential. 



Adjustable parameters 

We did not adjust parameters to try to achieve the best 
possible fit to the data. Rather, we chose values for pa- 
rameters as close as possible to the values in the real 
device while maintaining computational feasibility. 

Values for most of the important parameters for the 
DFT calculation are given in the manuscript. Other rel- 
evant parameters are the lithographic width of the quan- 
tum wire away from the QPC (300 nm), the voltage on 
the quantum wire electrodes (V c — — 0.08V) and the po- 
sition of the donor layer (20 nm below the surface). The 
details of our numerical approach are presented in a sep- 
arate publication [2ll |. 



• detector (in paper) 



Figure SI: a) The electrodes of the modified QPC used in the 
DFT calculation, b) The cross section through the center of 
the QPC (dashed line in panel a). The red dot 200 nm below 
the 2DEG is the position of the detector in the calculation 
presented in the paper. The derivative D DFT calculated from 
the potential in the 2DEG plane 300 nm from the center of the 
QPC (blue dot) is shown in the next figure. In the calculation 
presented in the paper the voltages on both gate electrodes 
were changed together (V^ pc ,L = VqpcFt), while in the present 
calculation Vqp C ,R was fixed and we only varied Kj PC ,l. 
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Can spin rearrangement be detected 
with our method? 

The additional dip near pinch-off in the local spin-density 
approximation (LSDA) calculation is related to forma- 
tion of a spin-1/2 magnetic moment in the QPC. Some- 
where between pinch-off and the transition to the first 
conductance plateau, a quasibound state forms at the 
center of the QPC just below the Fermi energy. The 
quasibound state can be occupied with either a spin-up 
or a spin-down electron, i.e. the QPC acts as a spin- 
1/2 magnetic moment. The formation of the magnetic 
moment is accompanied by a slight rearrangement of the 
electron density within the QPC, which shows up as an 
additional dip in D. 

As alluded to above, some of the present authors have 
performed an extensive study of the conditions under 



sitive measurements may be able to detect it. 



which such a magnetic moment forms in the QPC [21 1. 
We observed this feature only in the DFT calculation 
within the local spin-density approximation. In the ordi- 
nary local density approximation the two spin-densities 
are restricted to be the same and thus spin-polarization 
in the QPC is not possible, while in the Hartree approx- 
imation the exchange interaction, which is responsible 
for spin-polarization, is absent. We observed the forma- 
tion of the magnetic moment for various geometries of 
the QPC (rectangular shape of QPC electrodes or more 
adiabatic, triangular shaped electrodes) and for a wide 
range of QPC widths and 2DEG densities. The param- 
eter which affects the formation of the quasibound state 
in the most important way is the length of the QPC. For 
a substantial range of QPC lengths the situation is the 
same as presented in this paper. However, if the QPC is 
very short there is not enough space for the quasibound 
state to form in it. If the QPC is very long an antiferro- 
magnetically ordered spin chain forms instead of a single 
spin-1/2 magnetic moment. 

Our present experiment cannot resolve the extra dip 
described by the LSDA calculation, but future, more sen- 
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